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Abstract 

The multiparticle spacetime algebra (MSTA) is an extension of Dirac 
theory to a multiparticle setting, which was first studied by Doran, Gull 
and Lasenby. The geometric interpretation of this algebra, which it in- 
herits from its one-particle factors, possesses a number of physically com- 
pelling features, including simple derivations of the Pauli exclusion prin- 
ciple and other nonlocal effects in quantum physics. Of particular im- 
portance here is the fact that all the operations needed in the quantum 
(statistical) mechanics of spin 1/2 particles can be carried out in the 
"even subalgebra" of the MSTA. This enables us to "lift" existing results 
in quantum information theory regarding entanglement, decoherence and 
the quantum / classical transition to space-time. The full power of the 
MSTA and its geometric interpretation can then be used to obtain new 
insights into these foundational issues in quantum theory. A system of 
spin 1/2 particles located at fixed positions in space, and interacting with 
an external magnetic field and/or with one another via their intrinsic 
magnetic dipoles provides a simple paradigm for the study of these issues. 
This paradigm can further be easily realized and studied in the laboratory 
by nuclear magnetic resonance spectroscopy. 

1 The Physics of Quantum Information 

Information, to be useful, must be encoded in the state of a physical system. 
Thus, although information can have more than one meaning, it always has 
something to say about the state of the system it is encoded in. Conversely, the 
average information needed to specify the state of a system drawn at random 
from some known probability distribution is proportional to the entropy of the 
corresponding statistical mechanical "system". It follows that entropy can be 
understood as a measure of the system's information storage capacity. The 
physics of information is a fertile area of research which promises to become 
increasingly important as computers and nanotechnology approach the limits of 
what is physically possible |^ . 

In practice, information is usually binary encoded into an array of two-state 
systems, each of which can hold one bit of data, where the two "states" corre- 
spond to the minimum or maximum value of a continuous degree of freedom. 
Ordinarily such physical bits obey classical mechanics, but many examples of 
two-state quantum systems are known, for example adjacent pairs of energy 
levels in atoms, photon polarizations, or the magnetic dipole orientations of 
fermions. These quantum bits, or qubits as they are called, have a number of 
distinctive and nonintuitive properties P, E4, 03]. In particular, the number of 
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parameters needed to specify (the statistics of measurement outcomes on) the 
joint state of an array of n qubits is generaUy 2" — exponentiaUy larger than 
the n needed for an array of classical bits! These new degrees of freedom are due 
to the existence of non-separable or entangled states, which may be thought of 
as providing direct paths between pairs of states related by flipping more than 
one bit at a time. A further mysterious property of these quantum states stems 
from the fact that it is impossible (so far as is known!) to determine just where 
on the pathways between states the qubits are. This is because the act of mea- 
suring the qubits' states (however this may be done) always "collapses" them 
into one of their extremal states jl^ . 

A quantum computer is an array of distinguishable qubits that can be put 
into a known state, evolved under a sequence of precisely controlled interactions, 
and then measured. It has been shown that such a computer, if one could be 
built, would be able to solve certain problems asymptotically more rapidly than 
any classical device. Unfortunately, quantum systems are exceedingly difficult to 
isolate, control and monitor, so that at this time only simple prototype quantum 
computers have been operated in the laboratory. Although far from competitive 
with today's laptops, these prototypes are of great scientific interest. This is 
because quantum computers provide a paradigm for the study of a number of 
poorly understood issues in quantum mechanics, including why the particular 
classical world we inhabit is singled out from the myriads allowed by quantum 
mechanics. This is widely believed to be the result of decoherence: the decay 
of accessible information due to the entanglement generated by the interactions 
between the system with its environment ||, |l^ . The reason such issues are 
poorly understood, even though the microscopic laws of quantum mechanics are 
complete and exact, lies in our very limited ability to integrate these laws into 
precise solutions for large and complex quantum systems, or even to gain signif- 
icant insights into their long-term statistical behavior. The intrinsic complexity 
of quantum dynamics is in fact precisely what makes quantum computers so 
powerful to begin with! 

This paper will explore the utility of the multiparticle spacetime algebra 
(MSTA) description of qubit states, as introduced by Doran, Gull and Lasenby 
H 0, |g], for the purposes of understanding entanglement, decoherence and 
quantum complexity more generally. As always with applications of geometric 
algebra, our goal will be to discover simple geometric interpretations for oth- 
erwise incomprehensible algebraic facts. To keep our study concrete and our 
observations immediately amenable to experimental verification, we shall limit 
ourselves to the qubit interactions most often encountered in physical imple- 
mentations, namely the interaction between the magnetic dipoles of spin 1/2 
particles such as electrons, neutrons and certain atomic nuclei. These will be 
assumed to have fixed positions, so their spatial degrees of freedom can be ig- 
nored. Examples of such systems, often involving 10^° or more spins (qubits), 
are widely encountered in chemistry and condensed matter physics, and can 
readily be studied by various spectroscopies, most notably nuclear magnetic 
resonance (NMR) |2|. 
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2 The Multiparticle Space-Time Algebra 



It will be assumed in the following that the reader has at least a basic familiarity 
with geometric algebra, as found in e.g. ||^, ^, This brief account of the 
MSTA is intended mainly to introduce the notation of the paper, while at the 
same time providing a taste of how the MSTA applies to quantum information 
processing. More introductory accounts may be found in [Q, ^ 

The rt-particle MSTA is the geometric algebra Q{n, 3n) generated by n copies 
of Minkowski space-time TZ^'^. We let {7^ | ^ = 0, . . . ,3} denote a basis set of 
vector generators for this algebra, satisfying 

= v^.s'^'- (1) 

The superscripts refer to separate particle spaces (a, 6 = 1,... ,n) and the 
subscripts label spacetime vectors (/i, ~ 0, . . . ,3), while rj^i, is the standard 

Minkowski metric of signature (H ). The {7^} can be thought of as a basis 

for relativistic configuration space, and the MSTA is the geometric algebra of 
this space. 

The even subalgebra of each copy of the spacetime algebra is isomorphic to 
the algebra of Euclidean three-dimensional space ^/(3) The specific map 

depends on a choice of timelike vector, and the algebra then describes the rest 
space defined by that vector. It is convenient in most applications to identify 
this vector with 70, and we define 

f^k = 7fe7o • (2) 

Each set of spatial vectors (Ti,cr2,(T^ generates a three-dimensional geometric 
algebra. It is easily seen that the generators of different particle spaces com- 
mute, so that the algebra they generate is isomorphic to the Kronecker product 
t/(3) (g) ■ • ■ (8) 0(3). For this paper we will work almost entirely within this (non- 
relativistic) space, but it should be borne in mind throughout that all results 
naturally sit in a fully relativistic framework. 

To complete our definitions, we denote the pseudoscalar for each particle by 

= 7o7i72 73 = '^i'^2'^3- y-^) 

For the bivectors in each spatial algebra we make the abbreviation t^cr^ = «t^. 
The reverse operation in the MSTA is denoted with a tilde. This flips the sign 
of both vectors and bivectors in Q{'i), and so does not correspond to spatial 
reversion. 

There are many ways to represent multiparticle states within the MSTA. 
Here we are interested in an approach which directly generalizes that of Hestenes 
for single-particle non- relativistic states. We will represent states using multi- 
vectors constructed from products of the even subalgebras of each t/(3). That 
is, states are constructed from sums and products of the set {l,/cr|'}, where 
k = 1 ... 3 and a runs over all particle spaces. This algebra has real dimension 
4", which is reduced to the expected 2"+^ by enforcing a consistent representa- 
tion for the unit imaginary. This is ensured by right-multiplying all states with 
the correlator idempotent 

E ^ uriuri) i(l - «t>|) • ■ • i(l - /^jO . (4) 
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The correlator is said to be idempotent because it satisfies the projection relation 

= E. (5) 

In the case of two spins, for example, we have E — —ut^ut^E, so that any term 
right-multiplied by —ur^ur^ is projected by E back to the same element. The 
correspondence with the usual complex vector space representation is obtained 
by observing that (for e.g. two spins) every spinor may be written uniquely as 

i/) = - UT^xp2 - UT2ip3 + i^li^i^i) E , (6) 

where the i/)^ are "complex numbers" of the form + /3fcJ with a^, Pk real, 
and the role of the imaginary is played by the (non-simple) bivector 

J ^ UT^E = uriE = i(«T3i + «t|). (7) 

The complex generator J satisfies = — E, which ensures consistency with 
the standard formulation of quantum theory. As in the single-particle case, the 
complex structure is always represented by right-multiplication by J. While this 
approach may look strange at first, it does provide a number of new geometric 
insights into the nature of multiparticle Hilbert space H . 



3 Two Interacting Qubits 

As an application of the MSTA approach, we consider a simple model system 
of interacting qubits. This exhibits all of the complexity of multiparticle quan- 
tum systems generally, including the role of entanglement and the distinction 
between classical and quantum theories. Associated with the intrinsic angular 
momentum of any spin 1/2 particle is a magnetic dipole fi. In the far-field limit, 
the energy for the interaction between two such dipoles is given classically by 
the expression 

^o^£(^^^-3Hi^). ,8) 

where is the permeability of the vacuum, r is the radial vector between the 
dipoles, and r = \r\ (see Fig. |^). To obtain the quantum theory of this system 
(via "first quantization") we replace the magnetic moment by its operator 
equivalent, given by the component relation 

Afc = JSk- (9) 

Here 7 is the gyromagnetic ratio and Sk is the spin operator in the fc-th direction. 
For a spin 1/2 particle the spin operators are simply 

Sk ^flCTk, (10) 

where the at are the Pauli matrix operators. The classical energy of (^) gives 
rise to a quantum Hamiltonian containing 2 terms. The first involves fii - fj,2 , 
which is replaced by the operator 

EAlA^. = ^E-'^^^^' (11) 

fe=l k=l 
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Figure 1: One dipole in the field of another, the radial vector r between them, 
and the angle 9 it makes with the dipole (left); corresponding potential energy 
surface (assuming the dipoles are confined to the plane as depicted on the left) 
showing that there are two minima, two maxima, and four saddle points (right). 



where 71 and 72 are the gyromagnetic ratios of particles 1 and 2 respectively. 
This operator acts on a four- dimensional complex vector space. To form the 
MSTA equivalent of this we write 

3 3 
^ <Tfc <8) o-fe = - ^ i<Tft (8) iak , (12) 

where i is an (uninterpreted) imaginary unit. Each factor of i&k has an equiva- 
lent action in the MSTA given by left-multiplication with UTk in the correspond- 
ing particle's space. It follows that we can replace 

3 3 
^CTfeOo-felV) ^ -^ta-luT^ip . (13) 
fe=i fe=i 

We can already see, therefore, that the Hamiltonian is going to become a 4- 
vector in the MSTA. 

For the second term we let n = r/r, so that n is the unit vector parallel to 
the line between the dipoles' centers. We next form the operator for n^- n ^12- n, 
which is 

3 ^2 3 

^ nkTiefilfij = — - — ^ nkUiak^ut. (14) 
The MSTA equivalent of this is simply 

ft^lljl ^ 1 2 fl^lll2 12 

— 2^ nk ni OTfe /o> = — m m . (15) 

k,e=i 
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The role of the Hamiltonian in the MSTA is therefore assumed by the 4-vector 

^ - I (E - 3 . (16) 

k=l 

where d = /io^7i72/47rr'^ sec^^. This definition of Hq is chosen so that 
Schrodinger's equation takes the simple form 

~dtip - HD^pJ. (17) 

^From this one can immediately see a key feature of the MSTA approach, which 
is that all references to the tensor product have been removed. All one ever 
needs is the geometric product, which inherits all of the required properties 
from the relativistic definition of the MSTA. 



3.1 The Propagator 

In the more conventional, matrix-based approach to quantum theory the propa- 
gator in the current setup would be simply exp{—iHt/h). In finding the MSTA 
equivalent of this we appear to have a serious problem. In equation ( p^ the 
4-vector Hq acts from the left on tp, whereas J sits on the right. In fact, this 
is more of a notational issue than a foundational one. We simply define an 
operator j to denote right-multiplication by J when acting on a multiparticle 
state, namely 

jiP = iPJ. (18) 

We are then free to write j anywhere we chose within a multiplicative term, and 
to let it distribute over addition like a multiplicative operator. This notational 
device is extremely useful in practice, though occasionally one has to be careful 
in applying it. 

The propagator can be obtained in a number of ways. One is to write 

Ho = -n - - + -m}m.\ (19) 

where 

n = i (1 - uctIuct^ ~ ur^uri ~ uriurl) . (20) 

The multivector II constitutes a geometric representation of the particle inter- 
change operator, since UC^v^ = u^(^^n, and hence satisfies 11^ — 1. It follows 
that II commutes with <n^/n^, and the propagator can be written as 

expi-jHot) = ex.p{-j dt/4)exp{jUdt/2)exp{j m^m'^3dt/A). (21) 

All three exponentials in this expression commute. 

Alternatively, one can look for eigenstates of the Hamiltonian. At this point 
it is convenient to choose a coordinate system in which r is parallel to the z-axis, 
so that 

n = (T3 , Hd = -(2«T3/<t| - urlirrl - utIut^) ■ (22) 
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One can "diagonalize" the Hamiltonian operator in this coordinate frame by 
defining 

d / I!? '12/12 1 2\\ 

= ^ (2 «^>3 + J ("^3 - 1^1)) ■ (23) 
The eigenvalues and eigenspinors of H'^ are 

H'qE = —-E, H'qUT2E = duT2E , 

H'^uj^E = 0, if'D«r2W|£; = -'^ujIuj^E. (24) 

It follows that the propagator in the transformed basis can be constructed as 

+ e-J'^'^^t^^UT^EuT^uri) (25) 
= e-odt/2E - \{l + ei<'')uT^EuT^ , 
where we have used the results 

UT2UT2EUT2UT2 = E and UT2EUT2 = ut2EutI. (26) 
The propagator is easily transformed back into the original frame to give 
e-iJ/D t ^ £;_!(! + e^''**) {url E ur^ + «t| E ut^) 

+ i (1 - e^^') {uri E uri + «t| E ut^) (27) 
= e-^'*/2 E - ^{1 + e^''')uT^EuT^ + 5(1 " eJ''^*)«T2i£;/£r| , 

where we have made the same simplifications as before. 
3.2 Observables 

The simple form of the propagator hides a number of interesting properties, 
which emerge on forming the observables. Given an initial state i/>q the state 
at time t is simply 

^{t) = e-^'«-*Vo- (28) 

Assuming ||t/'oI| = 1, the spin bivector observable is 

S{t) = 2ip{t)J'ip{t), (29) 

which defines separate spin bivectors in spaces 1 and 2. This object turns out to 
have a number of remarkable features. If the two spins start out in a separable 
state with their spin vectors aligned with the z-axis, i.e. ipQ = E, then they do 
not evolve since Hq E = {—d/2)E, collapsing the propagator to a phase factor. 
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This result is in accord with the fact that this orientation is a classically stable 
equilibrium. 

A classically W77.stable equilibrium is obtained when the spins are aligned 
antiparallel along the z-axis, e.g. in the quantum state i/'q = /crj^i^. The time- 
dependent spinor in this case is given by 

2e~^-^''*«Tj£; = ((l + e^-^*)«Tj - {l - e^''')uri)E ^^^^ 
= (ttT2 - 10-2) E + (ttT2 + UT2) ( cos(di) E + sin(dt) j) , 

so that the state oscillates between UT2E and —UU2E with a period of iir/d. 
This spinor can be written in the canonical form (obtained by singular value 
decomposition pT|) 

e-j«Dt^^i£; ^ urle^"'^-"'^'^^''^ E , (31) 

from which it is easily seen that the spin bivector is 

S{t) = cos{dt) {utI - utI) . (32) 

Although they remain equal, the magnitudes of the two spins' bivectors can 
shrink to zero, showing that they have become maximally entangled. Thus we 
see that the state oscillates in and out of entanglement while swapping the signs 
of the spin bivectors every half-cycle. 

Now suppose that the spins start out with both their vectors parallel along 
the X-axis, i.e. with i/jq = ^(l + i^2){^ + '"'l)-^' which is a saddle point of 
the classical energy surface. Then our propagator gives us the time-dependent 
spinor 

e-J'-f^^'-^o = \({l + UTluT^)e-i'''/^ + {u^l + UTi)ei'''^E , (33) 
which in turn gives rise to the spin bivector observable 

S{t) = i K + >^2) J e^'"''' (1 + «T>|) - 

\ (1 4- «T>|) J e-^^'^'^ {^\ + u^i) (34) 
= cos(3dt/2) [utI + UT^) . 

It follows that the measure of entanglement in this case varies as the cosine 
oi •& = Sdt/2. Curiously, however, if the spins start out antiparallel, a similar 
calculation gives d = dt/2. Of course the dynamics are the same if the spins 
start out (anti)parallel along the y-axis. 

Finally, Fig. ^ shows plots of the trajectories of the spin vectors on (or near) 
the surface of a unit sphere, starting from an unentangled state with the first 
spin (light gray) along the x-axis and the second (dark gray) along the z-axis. 
The lengths of the spin vectors stayed very nearly at unity, implying that little 
entanglement was generated. The first spin executed a loop-de-loop up to the z- 
axis, returned to the —x-axis, and continued on in the general direction of the y- 
axis, while the second swooped down towards the y-axis, returned symmetrically 
to the z-axis, and then ended half-way between the x & z-axes. The complexity 
of the trajectories even in such a simple quantum system is impressive! 
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Figure 2: Plot of the spin vector trajectories under the dipolar Hamiltonian 
starting from an unentangled state with the first spin (light gray) along the x- 
axis and the second (dark gray) along the z-axis. The length of the vectors stayed 
near unity, implying that the state remained largely unentangled throughout its 
acrobatics. 

If we take the derivative of the spin bivector we get 

dtS{t) = 2{~HDipJ)J'ip + 2ipJ{-HDtpJT 

= 2{HdiP^ - tpi^Ho). (35) 

Now Hd is an MSTA 4- vector, while ipip is the sum of a 4-vector and a scalar, 
contains a 4-vector term as well as the scalar term. The commutator of two 
4-vectors gives rise to terms of grade two, and it is this that gives rise to the 
interesting dynamics. If the state is separable we can write 

2{i,i,)^ = qa'tq^ (36) 

where p and q are the single-particle spin vectors. In this case we have 

dtS{t) = Hoip^iq^ -ip^tq^Ho. (37) 

For our Hamiltonian the general relation for commutators of tensor products of 
three-dimensional bivectors, 

(A^B^) X (X'Y^) EE ^{A^B^){X^Y^) - 1{X'Y^){A^ B^) 

= A-X {B xY)^ + B-Y {Ax X)\ (38) 

enables us to reduce this to 

9tS{t) = ^(pAq)^-^(pAq)"+yO■3•q(a■3Ap)l + yO•3•p(o•3Ag)^ (39) 
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and we have recovered exactly the classical equation of motion for two dipolar- 
coupled spins. Of course, separability is not preserved in the case of quantum 
spins, due to differences in the higher derivatives. This allows the spin vectors 
to change their lengths (as the degree of entanglement varies), which cannot 
happen classically. 



4 Lagrangian Analysis 

Lagrangian methods play a central role in finding equations of motion for con- 
strained systems such as rigid bodies in classical mechanics. Via Noether's 
theorem, they also permit one to identify constants of the motion such as the 
total angular momentum, which provide insights into their long-term dynamics. 
Since quantum dynamics typically exhibit many more symmetries than classical, 
an extension of Lagrangian methods to multispin systems is highly desirable. 
This type of analysis is somewhat underutilized in traditional quantum treat- 
ments, as it sits somewhere between classical Lagrangian analysis and quantum 
operator techniques. 



4.1 Single-Particle Systems 

As a starting point, consider a single spin-1/2 particle interacting with an ap- 
plied magnetic field described by the bivector B. The Lagrangian for this is 



L = i^xl) UT3 - ba-^tl)) . (40) 

Substituting this into the Eulcr-Lagrange equations 

d_ fdlA _ dL_ 



yields 



(41) 



^{uTs'ip) = {■ipUTzy - -/{uTs'lpB) - j{BipUT3y . (42) 



This simplifies readily to Schrodinger's equation for the system: 

-ip ^ jBip. (43) 

This prototype can be used to build more realistic semi-classical models of spin, 
including relativistic effects 0] . A feature of this Lagrangian, which is typical for 
spin-1/2 systems, is that it is first order in xp. For such systems one frequently 
finds that L — for paths satisfying the equations of motion. 



4.2 Two-Particle Interactions 

Of greater relevance here is the two-particle equation (^^. This can be obtained 
from the Lagrangian 

L{ip,xp) = 2(ipJ'ip - HotpE'ip) . (44) 
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The first term here can be viewed as the kinetic energy, and the second term 
as the potential energy. The latter term couples entirely through the 4-vector 
part of ip E'ip. Some insight into the nature of this system can be obtained by 
parameterizing the wavefunction as 

ip = v^C'^{cos{e/2)E + sm{e/2)uT^uT^J)e°'-^. (45) 

Here v and C are single-particle spinors, 9 measures the entanglement and a is 
an overall phase factor. In total this parameterization has 10 degrees of freedom, 
so 2 must be redundant. One of these is in the separate magnitudes of v and ^, 
since only their product is involved in tp. The second redundant parameter is in 
the separate single-particle phases, since under the simultaneous transformation 

i,^i,e-'^'^, C^Ce'^'^, (46) 

we see that ip is unchanged. Despite this redundancy, the parameterization is 
extremely useful, as becomes clear when we write the kinetic term as 

2(^1} J ij}^ = COS{0){CUT3C i^v) + VUT3V (CO) ■ (47) 

This would reproduce the classical dynamics of two magnetic dipoles, were it 
not for the factor of cos{9). We also see that there is no derivative term for 6, so 
the Eulcr-Lagrange equation for the entanglement measure produces a simple 
algebraic equation. 

The potential term in the Lagrangian can similarly be written 

+ sm{6) (^Hd CV2(ttTiW| + uritcrDc'v^) ■ (48) 

Again, it is the entanglement factor that adds the quantum effects. If we set 
= in the Lagrangian then the equations of motion reduce to precisely those 
for a pair of classical dipoles. It is the presence of 9 in the Lagrangian, which 
is forced upon us by the nature of multiparticle Hilbert space, that makes the 
system truly quantum mechanical. 

4.3 Symmetries and Noether's theorem 

Symmetries of the Lagrangian give rise to conserved quantities via Noether's 
theorem. Let ip' = ■0'(i/',q:) be a differentiable transformation of the spinor ip 
controlled by a single scalar a, and satisfying ip'{ip,0) = ip. If we define the 
transformed Lagrangian as 

L'{^p,i,) = L(tA',i/;') , (49) 
we find (using the Euler-Lagrange equations) that 

dL' d / dtp' dL\ „ d / ^ ~ \ 

where cj) = dtp' / da\a=o . If L' is independent of a our transformation defines a 
symmetry of the Lagrangian, and gives rise to a conjugate conserved quantity. 
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As a simple example, take the invariance of L under changes of phase: 

V?' = V'e"-^ => ct) = ipJ . (51) 

It is easily seen that this is a symmetry of the Lagrangian, and that the conjugate 
conserved quantity is —2{ipip), telling us that the magnitude of is constant. 

Phase changes are something of an exception in that they involve operation 
on ip from the right. Operation from the left by a term of the form exp(aP) will 
always generate a symmetry provided P commutes with the Hamiltonian. The 
prime example in this case is the Hamiltonian Hq itself. The symmetry this 
generates corresponds is time translation, as defined by the exp(— jifoi) and 
the conjugate conserved quantity is the total internal energy. Another important 
symmetry generator is provided by ut^ + ut^ , which generates rotation of both 
spins about the z-axis at equal rates. (Recall that this axis is defined by the 
inter-dipole vector.) This bivector commutes with the Hamiltonian since 

Hd {ut^ + «t|) = - ^ {uT^ + «t|) = {uT^ + uri) Ho , (52) 

and the conjugate conserved quantity is 

{UT^ + «t|) • (2t/> JxP) = (iTTgi + «t2) • 5, (53) 

which gives the total angular momentum about the z-axis. 

Yet another symmetry generator is found by considering the operator for 
the magnitude of the total angular momentum, namely 

{utI + irr^f + (wrj + ur^f + {^l + «t|)' ^^^^ 
= 2 [uctIutI + UT2UT2 + utIutI - 3) = - 4 (1 + n) . 

The new generator here is jll, which generates a continuous version of particle 
interchange, and is a symmetry of L. The conjugate conserved quantity is 
(n'0'0). Rather confusingly, this is not the same as the total magnitude (SS), 
which is not conserved in the quantum case (though it is conserved classically 
where the spin vectors have fixed length). 

Together with phase invariance, we have now found a total of four physically 
relevant constants of motion for the two-spin dipolar Hamiltonian Hq . Are 
there any more? Clearly any linear combination of constants of motion is again 
a constant of motion, so this can only be answered in the sense of finding a 
complete basis for the subspace spanned by the constants of motion in the 16- 
dimensional space of all observables. The dimension of this subspace can be 
determined by considering the eigenstructure of Hq again (Eq. |2j). Clearly 
the observable which measures the amount of a state in any one eigenspinor 
is a constant of the motion, and more generally, each eigenvalue gives rise to 

constants of the motion where 5 is its degeneracy. Since Hq has one two- 
fold degenerate eigenvalue, the constants of motion must span a subspace of 
dimension 1 + 1 + 2^ = 6, and we are therefore just two short! It is quite easy 
to see, however, that the so-called (in NMR) double quantum coherences 

commute with Ho and generate two new independent symmetries. The physical 
interpretation of these symmetries is far from simple as they again involve the 
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geometry of the 4-vector {xp'ilj)^. Because the product of any two constants 
of motion is again a constant of motion, the constants of motion constitute a 
suhalgehra (of the even subalgebra) of the MSTA. This makes it possible to form 
large numbers of new constants in order to find those with the simplest physical 
or geometric interpretations. 

5 The Density Operator 

For a single-particle normalized pure state the quantum density operator is 

p = = \{i + T.LiVk(Tu), (56) 

where / is the identity operator on a two-dimensional Hilbcrt space, and 

Pk = (V' I CTfc I V') • (57) 
The MSTA equivalent of this is simply 

^Z+ip = i(l+p) , (58) 

where Z± = i(l ± cr^) are idempotents and p — cr^xp is known as the Block 
or spin polarization vector. 

The most straightforward extension of this to the multiparticle case replaces 
the idempotent by the product of those for all the n distinguishable spins: 
= ■ ■ ■ Z" . This has the added benefit of enabling one to absorb all 
the vectors from E into Z+ , thereby converting E into a pseudoscalar 
correlator, 

C = i(l-ti/.2)... (59) 

This commutes with the entire even MSTA subalgebra, and allows one to iden- 
tify all the pseudoscalars wherever they occur, with a single global imaginary 
as in conventional quantum mcchancs. This faithfully reproduces all standard 
results, but the idempotent Z^ ■ ■ ■ Z^ does not live in the even subalgebra, 
and hence mixes up the grades of entities which otherwise would have a simpler 
geometric meaning. In the following we propose for the first time a formulation 
of the density operator within the even subalgebra of the MSTA. 

The key is to observe that both Z^ ■ ■ ■ and E + J contain all 2" possible 
products of the (a = 1,... ,n). It follows that the even MSTA density 
operator 

p, EE 2''-^iP{E + J)i, (60) 

includes the same complete set of commuting observables as does the usual 
definition (modulo pseudoscalar factors). The normalization of 2"~^ is chosen 
so that the scalar part (p^) — 1. This is more natural when using geometric 
algebra than the usual factor of 1/2, which ensures that the trace of the identity 
state 7/2" is unity. The density operator of a mixed state is simply the statistical 

average of these observables as usual: = ip{E -\- J)xp |jl^, pT[ |. 

There is a problem with the even MSTA version, however, which is that 
xp Jxp is anti-Hermitian (under reversion) , whereas the usual density operator 
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is Hcrmitian. In addition, up to now wc have only applied the propagator 
to spinors, where the complex structure is given by right-multiplication with 
J. Such a representation is not appropriate for density operators, because to 
apply it we would have to decompose our density operator into a sum over an 
ensemble of known pure states, apply the propagator to each state's spinor, and 
then rebuild the density operator from the transformed ensemble. This is clearly 
undesirable, so we should re-think the role of the imaginary in the propagator 
for the density operator. 

To do this, we must first understand the role of i in the propagator, where 
it appears multiplying the (Hermitian) Hamiltonian H . The terms in iH in- 
volving odd numbers of particles have an immediate counterpart in the MSTA 
as products of odd numbers of bivectors. These exponentiate straightforwardly 
in the MSTA, and here the role of the imaginary as a pseudoscalar is clear. It 
is the terms in iH involving even numbers of particles which are the problem. 
When these are converted to products of even numbers of bivectors, a single 
factor of i is left over. To see its effect we return to a pure state and consider 

{jtp)ij^ = tpJ'ip and {jil))J'^ = —tpij}. (61) 

Thus we see that, applied to density operators (or any other observables) , the 
imaginary unit interchanges even and odd terms (in their number of bivectors) , 
in addition to squaring to —1. 

It follows that the left-over factor of i in the even terms converts them to odd 
terms, thereby ensuring that the Hamiltonian generates a compact group while 
at the same time labeling these terms to keep them separate from the original 
odd terms. This is quite different from other occurrences of the imaginary! 
A simple way to represent this in the even MSTA is to introduce a "formal" 
imaginary unit j, and to redefine the even MSTA density operator as 

P = P+ - jp- (62) 

where = (pjg -f- (pj^ H and p = (pj^ + (pj^ H . Then, provided 

any residual factors of the even MSTA operator j are replaced by j, everything 
works simply. In writing p in this way we have also recovered a more standard, 
Hermitian representation, which permits the mean expectation values of any 
observable o = o+ + jo- to be computed by essentially the usual formula: 

{op)o = {{o+ + jo-){p^ - j p_))q = {o+p^ + 0-p_)o . (63) 

The geometric interpretation of this new imaginary operator j will be further 
considered on a subsequent occasion. 

5.1 An Example of Information Dynamics 

We are now ready to see how things look to one spin when no information about 
the other(s) is available. Clearly this depends on the (unknown!) configuration 
of the other spin(s), and we expect that the "generic" case will result in a 
very complicated trajectory, which may also be incoherent in the sense that the 
length of the spin vector is not preserved. Somewhat surprisingly, the situation is 
simplified substantially by assuming that the spin of interest has an environment 
consisting of a great many others, whose dynamics are so complex that their 
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spin vectors can be treated as completely random. Under these circumstances, 
neither the exact starting state of the environment, nor any subsequent change 
in it, will change the way our chosen spin (henceforth numbered 1) sees it. As 
a result of this assumption and the linearity of quantum dynamics, we need 
only figure out how its basis states ttr/ , UT2 , ur^ evolve in order to determine its 
evolution in general. 

To see how the kt^ evolve, we return to the propagator of equation (pT|), 
which we now write in the form 

exp(-j\ffDi) = ex.p{-juT^UT^Sdt/4:)eyip{-jUdt/2)ex.p{jdt/A) . (64) 

The phase term has no effect, and for the interchange term we find that 

e'-n^*/^ = cos\dt/2) ur^ + sm\dt/2) ur'^ 

- j sm{dt) n X tcr^ . (65) 

We then need to apply the term in exp{j ur^ ter^ 3dt / A) . This commutes with 
n, so we only need to transform the separate bivectors. For the case of fc = 3 
this final term has no effect, and we have 

e-^"°* U7^e3"°* = i(l + cos(dt))ttT3i + i(l-cos(dt))ttT| 

— j sin(di) i(«Tj^«T| - urluriy (66) 

For k — 1 , 2 we need the result 

g-j «T>3^0/2^_igj«T>3=^0/2 ^ cos((/))«Tj - j sin(0)/£ritcr|. (67) 
Combining this with the preceding we see that 

g-J «D t ^1 giHot ^ (cos(d<) UT^ - sm{dt) j ut^ut^) cos{dt/2) 

~ (cos(di) j utIutI + sin(dt) tcrl) sin(dt/2) , (68) 

with a similar result holding for ut\. 

As a simple model, suppose that our initial state is described by a known 
state of particle 1, encoded in the density operator p(0) = (1 +p(0))/2, and 
with the state of the second spin taken to be totally random. To evolve the 
density matrix of particle 1 we write <p(0) = ^Ylik Vk>t^\ and evolve each term as 
above. We then project p(t) back into the first particle space (by throwing out 
any terms involving other spins) and reform the polarization vector. We then 
find that: 

pit) = cos(dt) cos(dt/2) (pi cTi + p2 (72) + cos^(rfV2)P3cr3 • (69) 

Of course, the same result is obtained if everything is rotated to some other 
orientation, for example if the other spin is along x and p is replaced by the 
rotated spin vector. 

The von Neumann entropy of the spin, which measures how much informa- 
tion about its state has been lost to the environment, is given (in bits) by 

5vN = -i(l + ||p(t)||)log2((l + ||p(t)||)/2) 

-i(l-||p(t)||)log,((l-||p(t)||)/2) . (70) 

The signed lengths of the spin polarization vectors and the corresponding von 
Neumann entropies are plotted in Fig. ^ for a second spin in a random state 
and displaced from the first in directions parallel and perpendicular to its po- 
larization vector. 
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Figure 3: Dipolar "tides" (in units of 1/d) exerted on a spin's polarization vector 
by its interactions with a random ensemble of other spins all located in a circular 
"orbit" about its equator (light gray line), or at a fixed distance above its poles 
(dark gray line), along with the corresponding von Neumann entropies (dashed 
lines). 



5.2 Towards Quantum Complexity and Decoherence 

Our applications of geometric algebra to the problem of two dipole-couplcd spins 
has illustrated how quantum interactions allow qubits to exchange information 
(states) , in the process passing through an entangled state in which no informa- 
tion on either is directly accessible. We have also seen how information on the 
state of a qubit can be lost through its interactions with an environment, al- 
though the Poncaire recurrence times for the simple one-qubit environments we 
have considered are short enough to allow us to also see the underlying periodic 
behavior. This constitutes the simplest case of an old and venerable problem in 
solid-state NMR spectroscopy, which is to predict the spectral line-shape (decay 
envelope) of a macroscopic system of dipole-coupled spins such as calcium fluo- 
ride 19 . Here an exact treatment is out of the question, since the analytic 



evaluation of even the three-spin propagator is a reasonably challenging (though 
solvable) problem. Nevertheless, we believe the insights provided by geometric 
algebra, particularly into the constants of the motion, offer the hope of further 
progress. 

More generally, quantum mechanical explanations of apparently irreversible 
processes are becoming increasingly important in many fields [ p2[ , and further 
have the potential to finally clarify in just what sense classical mechanics can be 
regarded as a limiting case of quantum mechanics jo) . These problems typically 
involve the spatial degrees of freedom, and hence are truly infinite dimensional. 
Most of what is known about decoherence in such systems therefore comes from 
the analysis of simple and highly tractable environmental models such as a 
bath of spins or harmonic oscillators. Because of its potential to provide global 
insights into the full spacetime structure of such models, we expect that the 
MSTA will also come to play an enabling role in our understanding of the 
quantum mechanical mechanisms operative in decoherence. 
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